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CIRCULAR  RING  E-PLANE  EEEMENT  PATTERN  INTEGRAE  REPESENTATION 
AND  ITS  ASYMPTOTIC  EVALUATION 

Aly  E  Eathy,  A.  Hessel 

The  University  of  Tennessee,  ECE  Dept.,  Knoxville,  Tennessee  37996 

Abstract: 

Asymptotic  evaluation  of  the  far  field  for  large  kro  results  in  an  analytical  expressions, 
which  provides  a  significant  physical  insight  into  the  physical  process-taking  place  in 
element  pattern  formation. 


I-  Representation  of  a  Single  Element  Excitation 


Eig.  (1)  schematically 
illustrates  a  circular  array 
of  N  open  ended  parallel 
plate-waveguides,  in  a 
perfectly  conducting 
cylinder  coated  with  a 
dielectric  cover.  The 
waveguides  are  assumed 
to  propagate  solely  the 
TEM  mode.  Only  the  i  = 
0  element  is  fed,  the  others 
are  match  terminated. 


The  integral  representation,  regards  an  infinitely  extended  space  -oo  <  ^  <  oo,  in  which 
the  271  periodicity  is  enforced  by  image  elements  and  periodic  image  sources,  the  latter 
placed  an  angular  distance  27r  apart.  In  this  infinite  angular  space,  the  actual,  physical 
single  element  excitation  is  a  superposition  of  the  excitation  of  individual  sources  (s=0, 
±1,...).  Each  source  is  represented  in  terms  of  a  band  spectrum  of  angular  wave 
numbers  v. 


J-  =  0 


Figure  1  Cross  section  of  a  ring  array 


Denoting  the  angular  electric  field  component  by  E,|,(r,{zi;  v)  that  would  result  from  such 
an  excitation,  i.e.  the  field  due  the  excitation  of  the_7=  o  element.  In  the  feed  waveguide, 
the  modal  voltages  and  currents  are  defined  in  the  usual  manner.  The  field  at  an 
arbitrary  r  in  the  external  radial  waveguide,  can  be  expressed  in  terms  of  the  field  at  r  = 
ri,  as 


EAr,(l),  v)  = 


y: 


T(v.v,)  y 


(1) 


where  Vni=v+27r/pm,  and  P=27i/N.  Expression  (1)  relates  the  external  field  resulting 
from  phase  sequence  excitation,  to  the  incident  wave  voltage  Vine,  via  the  transmission 
coefficient  T(v,Vo)  and  the  voltage  ratio  V‘'(r^,vA/ V‘(r^,vA,  both  of  which  are 
dependent  on  the  form  of  aperture  electric  field  but  not  on  its  on  its  complex  amplitude. 
Y  o  is  the  admittance  of  the  fundamental  guide  mode,  and  Ga  is  the  real  part  of  the  active 
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admittance  Ya. 


II-  Element  Pattern 

Expressions  for  the  field  resulting  from  a  single  element  exeitation  are  obtained  with  the 
help  of  (1),  and  the  use  of  periodieity  in  the  v  of  the  Floquet  fields,  where 

71  CO  O 

^ ^  J ^  /*co 

=  and  in  view  of  the  identity  ^  f  {v  +  —m)dv  =  yf  {v)dv , 

jd  m=-<x>  H 

the  infinite  sum  of  the  finite  integrals  in  (1)  may  be  replaeed  by  a  single  infinite  integral 
[1,2,3],  i.e. 


2w,\G^{Voy^ 


nv,vj 


(2) 


As  k^r  ^  CO ,  even  though  p  <  oo ,  it  is  permissible  in  (2)  to  use  direetly  the  large 


argument,  first  order  asymptotie  expression  for  the  derivative  of  {kj)  .  Therefore, 
the  far  field  eorresponding  to  (2)  beeomes 


HV  (Vi) 


dv  (3) 


III-  Transmission  Coefficient  Tlv,  p,) 

III-l  Active  Admittance: 

When  the  form  of  the  aperture  field  is  a  slowly  varying  funetion  of  v,  the  behavior  of  the 
transmission  eoeffieient  determines  the  element  pattern  and  henee,  is  of  a  very  speeial 
importance.  Where  the  transmission  eoeffieient  is  funetion  of  the  aetive  refieetion 
eoeffieient  ra(v)  given  by: 

r(p,pj  =  i+r,(p,pj,  r,(p,pj  =  (7,(pj-7,(p))/(7:(pj+7,(p))  (4) 


From  (4)  it  is  seen  that  the  aetive  admittanee  Ya(v)  uniquely  speeifies  T(v,Vo)  for  a  given 
mateh  point  Vq  .  An  approximate  expression  for  Ya(v)  ean  be  obtained  as  in  planar  arrays 
from  the  requirement  of  eontinuity  of  complex  power  at  the  aperture  in  a  unit  eell  [1,2]: 


m=-co 


V 


p=i 


I  KM  K 


(5) 


where  Y*^  and  Y'  are  respeetively  the  modal  admittanee  of  the  radial  and  TFM  guides. 
Sinee  and  the  voltage  ratios  in  (5)  depend  on  the  form  of  the  aperture  eleetrie  field,  Eap, 
Ya(v)  may  be  ealeulated  only  provided  this  form  is  known.  For  the  ease  of  a  single  mode 
aperture  mode  will  be  treated  as  an  example.  The  first  term  in  the  right  hand  side  of  (5) 
represents  the  eontribution  of  the  exterior  region  to  the  input  admittanee  and  is  strongly 
v-dependent  beeause  of  the  form  of  the  dominant  radial  mode  eharaeteristie  admittanee 
Y‘^(r  ,v)  at  r  =  ro  whieh  is  is  given  by  [3] 


r  (r„,  1/)  =  V)  v)  =  -^  r'  (r,  k)  = 

Mo 


-J 


Hlikr^ 

MikrJ 


(6) 
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where,  k  =  k„Js,  and  the  air-dielectric  interface  reflection  coefficient  F  is 


r„  =  y^  ,7) 


yf^Y^,{kr,)  +  Ylik^r^) 


K^’ikx) 


K'^kx) 


III-2  Zeros  of  T  (v,  vj  and  Surface  Wave  Resonance  of  an  Unperforated  Cylindrical 
Structure: 


1.1 


•o.oe 


< 


-0.12 


-O.K 


T — I — I — I — I — I — I — 1 — r 


1.2 

—r~ 


l.<  R4vA  r  ) 


T — I — m — I — t — I — I — i 


T(v,Vo)  exhibits  zeros 
at  the  poles  of  Ya(v). 

These  occur  at  the 
value  of  V  for  which 
Y\r^,v,n)  =  ^-  These 
surface  wave  roots  for 
imperforated  structure 
occur  for  complex 
values  of  v  [1,2].  By 
continuity  argument 
with  the  planar  case, 
this  type  of  resonance 
is  expected  to  occur  Figure  2  Root  Locations  as  function  sheet  thickness 

for  a  value  of  Real(v)  for  which  the  air  region  is  "cut-off  and  the  dielectric  propagating. 
In  such  case,  kri>  Re(v)  >kori.  As  was  pointed  out,  these  poles  do  not  lie  on  the  real  v- 
axis,  as  in  the  planar  case,  since  an  attenuated,  angularly  propagating  surface  wave  is  not 
possible  in  the  cylindrical  case. 


•  t/2  .  0.2 

*  tA  •  o.is 
■  tA  •  0.1 


It  manifests  itself  in  a  resonant  behavior  of  the  imaginary  part  of  Ya(v),  where  the  latter 
experiences  a  sharp  rise  and  a  sign  reversal  (but  does  not  become  infinite  like  in  the 
planar  case)  in  the  predicted  range  for  a  value  of  Re(v/kori),  almost  equal  to  that  for  the 
corresponding  planar  case.  The  location  of  this  resonance  depends  on  the  dielectric 
sleeve  thickness.  It  is  seen  that  if  (t/k)  is  small  the  resonance  moves  into  transition 
region  where  (v=kori).  While  thicker  dielectric  covers,  locate  this  surface  wave  in  a 
region  where  Debye's  approximation  is  valid  for  the  different  cylindrical  functions.  The 
modal  admittance  reduces  to: 


ye'  ^  ^  -jcosw,^_ 

cos  ^  cos  j  cos 

where  j  =  -j  cot(^^,_  -  ^ )  and 


(8) 


4,  =kr{(x&w^  -(^-n4.)sinwj,.),sinw^  =  vl  kr  (9) 

Hence,  in  the  Debye  region  the  dispersion  relation  (8,  9)  for  the  surface  wave  on  smooth 
conducting  cylinder  covered  with  a  dielectric  sleeve  reduces  to 
cos  -  cos  =0  (10) 
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Fig.  2  shows  typical  root 
location  of  (10)  vs.  sleeve 
thickness  in  the  complex 
v-plane.  In  the  limit 
00  reduces  to 
relations  previously 

derived  in  the  planar  case, 
i.e.  to  the  lowest  TM 
surface  wave  on  a 
grounded  dieleetrie  slab  as 
was  pointed  out  earlier. 


!m(  v> 


N 


0.5 

too 


III-3  Poles  of  T(v,Vn): 


Figure  3  Pole  Singularities  of  T(v) 


From  (4),  the  poles  of  T(v,Vo)  (see  Fig.  3)  are  the  natural  resonance  of  the  mateh- 
terminated  array  structure,  and  correspond  to  the  periodic  solutions  of  the  transverse 
resonance  equation  where 

Y:{vJ  +  YSvO  =  ^,vi  =  y:+rn27,l  p,{m  =  Q,±\,...Up  =  ^,±h-)  (H)- 


The  pole  location  of  the  match  terminated  cylindrical  array  structure  is  a  perturbation  on 
that  of  the  short-cireuited  strueture.  The  dispersion  relation  of  the  later  corresponds  to 
the  poles  of  the  m  =  0  term  in  Ya(v),  i.e.  Equation  (11)  has  two  relevant  sets  of 
significant  poles.  One  arises  by  perturbation  due  to  finite  eurvature  of  the  planar  TM 
guided  waves  on  a  planar  match  terminated  array  structure.  The  other,  the  familiar 
creeping  waves  whieh  are  a  perturbation  of  those  on  a  full  dielectric  cylinder,  by 
bringing  the  perforated  conducting  cylindrical  surface  into  proximity  of  the  air-dielectric 
interface.  The  loeation  of  the  latter  is  near  v  =  kori. 


In  the  planar  case,  the  wave  numbers  of  the  corresponding  transverse  resonance  equation 
were  found  to  correspond  to  leaky  waves  and  the  fundamental  slow  harmonic  wave 
number  was  found  to  be  near  the  surface  wave  number  of  the  short-eireuited  array 
structure  shows  various  pole  locations  for  different  dielectric  thickness.  In  the 
cylindrical  case,  a  similar  situation  exists  for  the  first  set  of  significant  poles  where,  for 
the  set  m  =  0  of  (12)  the  relevant  surface  wave  pole  is  in  a  region  of  validity  of  Debye 
approximation.  In  this  case,  to  first  order,  the  modal  admittanees  are  identical  to  the 
corresponding  planar  expression.  Since,  the  creeping  waves  pole  locations  are  in  the 
transition  region,  i.e.  of  ,  v=kori  and  different  modal  admittance  representation 

must  be  used 


r'  = 


cosw, 


kr„ 


-^^Ar{T)Z  +  je 


-JYL  Ai\T)  +  jeA  {^Y^ZAi(T) 

cosw,^_  2 


(12) 


.Iti 


where  v  =  k/^  +  (■^^)'^^  re  ^ 


V  _  1  -  e 


sinw^^  =  — ,  Z  = 


kr. 


l  +  e 


-j2(?kn-?krA 

-j2(g,  -Ctr)  ^  J  kr,  ~‘^kr) 


April  19-23,  2004  -  Syracuse,  NY  ©  2004  ACES 


20th  Annual  Review  of  Progress  in  Applied  Computational  Electromagnetics 


where  Ai  is  Airy  function  and  ^kr  is  given  by 
=  kr{cosw^^  -  . 


IV-  Asymptotic  Evaluation  of  Far  Fields 

The  integral  representation  may  be  evaluated  asymptotically  for  large  values  of  kri 
yielding  simple 
analytical 

expressions.  The 
asymptotic  evaluation 
of  (3)  proceeds  in  a 
manner  similar  to 
smooth  large  cylinder 
or  uncovered 

cylindrical  array  [1]. 

Different 

representations  are 
obtained  for  the  lit 
region,  the  shadow 
region  and  the 
shadow  boundary 

(see  Fig.  1).  By  virtue  of  the  Physical  space  -n  <  ^  <  ;r ,  is  in  the  shadow  region  of  the 
image  sources  5  0  .  This  shadow  region  is  in  general  considerably  deeper  than  that  for 

s=0  source  and  consequently  the  contribution  from  s^o  are  neglected  except  near 
(p  =  ±n ;  (details  follow  later).  Flence,  the  lit  region  and  the  shadow  boundary  transition 
region  are  pertinent  only  to  the  s  =0  source. 


VrrCF.  IUi.LrfiJ4^i;q>-j 


^■4.  ii.  ■■■ 


A  Lit  Region 

The  basic  integral  in  (12)  to  be  evaluated  is  of  the  form 

¥{,!>)  =  r  T{y,  pjr  (q,  p)  / 

J— CO 


(13) 


under  the  assumptions  that  the  saddle  point  of  this  integrand  lies  in  region  where  Debye 
approximation  [1]  of  can  be  used  (this  assumption  will  be  justified  later),  i.e. 

\v^- k^r\>  0{k^ry^'^  hence,  the  integral  (13)  becomes 

V\r,,v) 


o'l 

nkf,  j 


J  T{v,v^) 

J— 00 


^dv 


cosw, 


(14) 


Kr, 


where  sinw^^  =  v !  kr^  and  ^^/'(p)  =  (^zi->v)sinw-cos>v 


For  a  given  value  of  (j),  the  saddle  point  Vs  or  Ws  are  located  at  /  t/p  =  o ,  i.e.  at 

<l)  =  w^  orsm(/)  =  smw^  =  vJ  kj^  (15) 

Since  the  range  of  Ws  is  limited  to  -;r/ 2  <  w  <  ;r/ 2 therefore,  has  a  dy/ !  dv  =  o .  This 
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condition  assures  that  Vg  <  kori,  and  justifies  the  assumption  of  validity  of  Debye 
approximation  for  This  solution  is  valid  only  in  the  lit  region.  Thus  the 

steepest-deseent  path  (S.D.P)  through  the  saddle  point  is  defined  by: 
Re  p)  =  Re  p  J  =  -  eos  (16) 


The  deformability  of  the  original  eontour  into  the  S.D.P.  in  the  v-plane  may  be 
aseertained  by  investigation  of  the  behavior  of  the  integrand  as  |p|  ^  oo .  By  Cauehy's 
theorem  one  has  (see  Fig.  4) 

=  {}r/p  +  2;z7-^[Re5(-pp-Re5(p;)]  (17) 

In  eontrast  to  smooth  eylindrieal  struetures,  the  appearanee  of  pole  eontributions  in  the  lit 
region  is  due  to  the  angular  periodieity  of  the  array.  The  evaluation  of  the  integral  in  the 
lit  region  is  eomplieated  by  the  faet  that  the  transmission  eoeffieient  may  have  poles  near 
or  on  the  S.D.P.  (Fig.  4).  In  the  development  to  follow,  the  first  order  saddle  point 
eontribution,  the  pole  residues,  and  the  transition  pattern  funetion  due  to  the  presenee  of 
a  pole  near  the  saddle  point  are  obtained. 


1:  First  Order  Saddle  Point  Contribution 

Upon  substitution  of 

F(Q)  =  J  f  (x)e  dx  =  yJlTr/Cl^q  (x)|/(x^)e 

Q^oo  for  q"  (x)  0 

the  first  order  S.D.P  eontribution  (direet  ray)  to  (13)  beeomes 


->4 


2  yUJpJ  F(r„,k„riSin^zi)  ) 


Tikr 


-j{kf,r-7r/4-k^rY  cos^) 


(18) 


(19) 


This  eontribution  is  observed  solely  in  the  lit  region,  appears  to  originate  from  the  center 
of  the  eylinder  and  is  referred  as  the  "direet  eontribution"  or  "the  spaee  wave". 

2,  Pole  Contributions:  From  (41),  the  residue  at  a  typieal  fast  wave  pole  of  vl  of  T(v) 


RQsivl)  =  e^Vfri,vl)/V‘ir^,vl)^7rk^rJ2cosw^,^a^ 
while  the  Residue  ap  is  given  by  =  hm_^^^,  (p  -  ppr(p) 


(20) 

(21) 


In  (20),  the  pattern  (j)  dependenee  enters  exelusively  through  the  exponential;  the 
remainder  of  the  expression  depends  only  on  the  pole  loeation  and  represents  the 
exeitation  eoeffieient  for  surfaee  waves  or  ereeping  waves.  The  poles  of  T(v)  are  ±p^ 

and  for  q=  0  represent  the  harmonies  of  the  fundamental  q  =  0  set.  With  the  help  of  (26) 
and  (47),  a  typieal  fast  wave  pole  eontribution  to  the  radiated  eleetrie  field  beeomes 


a^e 


-j(Kr+KhV(<)) 


(22) 


As  seen  from  Fig.  9,  only  the  q  <0  poles  are  eaptured  in  the  lit  region.  For  element 
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spacing  V2<a/A,<l,  the  poles  corresponding  to  q  <  -1  (2)  lie  in  the  region  |Re(v/kori)|>l 
where  exponentially  increasing  with  |v|.  Henee,  the  dominant  pole  eontributions  in  the  lit 
region  [1,2]  are  those  eorresponding  to  q  =  -1,  for  whieh  |Re(v/kori)|<l,  and  for  Wi  is 
small.  The  exponent  in  (49)  is  found  to  be 

-JKw(- )  =  -J'(^ - |)| I  +  +  Ah  cos w,  +\Vi\{<p-\w^\)  +  j-  (23) 

where  R  =  +  7  ,  w  =  w^+  jw.  and  for  (-p^'),  p^  >  0,  v.  >0,w^>  0,  w.  >  0)  . 

A  similar  phase  expression  for  those  poles  loeated  at  y=  and  is  given  by 

-jkjiWiy~p  )  =  ji</>  +  k.|)k.|+  cosw^  -  |u..|(^  +  |wj)  + 


Figure  5  Ray  Contributions  in  the  lit  Region 


where  Vr  <o,Vi,<o,  Wr 
<o,Wi  <0  for  poles. 

The  first  term  (23) 
expresses  the  faet  that 
with  inereasing  (J),  for 
(jp-O  the  effeet  of  this  set 
is  more  pronouneed 
(undergoes  less 

attenuation),  and  has 
maximum  eontribution 


a  relevant  pole,  i.e.,  {(j)  =  Wr  (whieh  is  the  pole  transition  region).  On  the  other  hand,  (24) 
shows  that,  this  set  has  less  effeet,  more  attenuation,  as  ^  is  inereased  for  ^zi>o.  For  the 
ease  of  lower  half  plane  poles  (w,<o,  Wr  <o,  Vi.<o,  Vr<o)  the  real  part  of  (24)  indieates 
that  a  typieal  ray  travels  eloekwise  through  an  angle  {(j)  +  |wr|)|  along  the  surfaee  and 
deeays  at  a  rate  of  |vi  |>  nepers/radians.  The  first  term  in  the  imaginary  part  of  (24) 
represents  the  phase  ehange  as  the  phase  ehange  as  the  ray  travels  from  the  souree  along 
the  surfaee;  its  assoeiated  phase  veloeity  is  sin  (wr)<0,  and  the  ray  thus  represents  a 
baekward  fast  wave.  The  seeond  term  in  the  imaginary  part  of  (24)  is  the  phase 
refereneed  to  the  eylinder  eenter,  assoeiated  with  the  free  spaee  path  of  the  ray  and 
implies  that,  irrespeetively  of  (j),  the  fast  set  emerges  from  the  surfaee  at  a  eonstant  angle 
Wr  with  the  loeal  normal  (the  positive  angles  are  eloekwise),  as  shown  in  Fig.  5. 


Rays  due  to  upper  half  plane  poles  follow  a  path  whieh  is  the  mirror  image  of  that  for  the 
lower  half  plane  poles,  and  (23),  (19)  show  that,  this  set  of  poles  has  a  phase  differenee 
with  that  of  the  direet  ray  given  by  lS.(/)Kkj-^{(X)S(/)-CQ&w^)  +  vJkj-^{(/)-w^),  where  Wr 
>o,Wi>o  and  very  small,  v,->o.  Fast  Wave  Pole  Transition  When  the  integral  along  the 
steepest  deseent  path  is  approximated  by  the  saddle  point  eontribution,  this  overall  result 
experienees  a  diseontinuity  as  a  funetion  of  (j)  every  time  the  path  erosses  a  pole,  the 
diseontinuity  being  proportional  to  the  pole  residue.  For  poles  deep  in  the  eomplex  plane, 
this  residue  is  exponentially  small  and  the  transition  effeet  may  be  negleeted.  For  poles 
elose  to  the  real  axis,  their  proximity  to  the  S.D.P.  may  be  taken  into  aeeount  and  a 
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continuous  transition  field.  The  overall  result  of  the  asymptotic  evaluation  of  the  integral 
(13)  is  given  by  the  sum  of  the  saddle  point  eontribution  (19),  the  residue  contribution 
from  the  fully  captured  poles  (22)  and  the  contribution  from  the  transition  poles.  The 
latter  harmonic  q  =  -1  set,  passes  through  this  transition,  where  S.D.P.  erosses  the  surface 
wave  root  ( )  first,  then  (  )  for  p>  0  and  the  eontribution  is  given  by 


Fm 


transition  pole 


=  e 


cos^T— ) 


) /  r ( vl )a^ ^Tdcjylcoswlx 


sin%=7^;^=-l  and  Xp  a/Vi 

Vi 

Q(y)  is  the  transition  function  given  by  the  complementary  error  function 
Q{y)  =  J  e  dx  ,  and  +  coscp 


(25) 


(26) 


B.  Shadow  Region  7t:/2<  I  <71 

The  integral  {nl2<(/Kn)  can  be  expressed  as  a  residue  series  at  the  poles  of  T(v).  For 
element  spacing  V2<a/A,<l,  the  eontributions  in  the  shadow  region  eome  mainly  from  the 
sources  s  =  0,  and  s=  ±1.  The  eontributions  of  all  other  sourees  are  highly  attenuated,  as 
they  have  to  travel  long  angular  distances  to  reaeh  the  physieal  spaee  -tkc/kti.  For 
example,  forF(^)  =  2;7g''^Re5(pp +  2;5r^Re5(-pp  .  Where  for  s=0,  the  path  of 

integration  is  elosed  at  the  lower  half-plane,  and  for  s  =  1,  it  is  closed  at  the  upper  half 
plane,  as  shown  in  Fig.  6.  A  elear  pieture  is  eaptured  by  studying  the  phase  for  the  slow 
wave  sets  =  0,  q  =  0, 

where;-7ViKv'°)  «  -\v\{(p  +  \w\-K)-jk^r^{cosw^+\v\  /  kxy.{(p-n+\w^\)-jKl  4  (27) 

but,  n !  2  and  >Vr<0,  Vr<0  .  So  for 

those  slow  waves  in  the  lower  half-plane 
one  has  Vi<0,  Vi>o,  Consequently 
-jKrM.^p)  ~  H  2)  -j  y.{(p-nl  2)  -jn!  ^ 

The  first  term  of  (27)  gives  an  attenuation 
as  if  it  were  traveling  an  angular  distance 
(XP -^\w n  ^  (p  -  n  !  2')  with  an 

attenuation  factor  |vi|.  A  corresponding 
relation  is  found  in  the 

phase,  i.e,  the  seeond  term  in  (27),  which 
shows  a  forward  wave  travels  a  positive 
angle  {(p +  n  (p  -  n  !  2')  from  the  excited  element  with  a  positive  phase  veloeity, 

as  shown  in  Fig,  11b,  while  the  third  term  Okoi’i  cos  Wr  ~  0),  so  when  the  phase  is 
eombined  with  (-jkor),  indicates  that  the  wave  departs  from  the  surface  at  a  point  B  (to 
whieh  it  travels  on  the  surface  an  angular  distanee  »  (j>-’Kl2  and  radiates  to  the 
observation  point  at  the  angle  (7i-|Wr|  «  7r/2)  with  respeet  to  loeal  normal).  For  fixed  p, 
the  higher  terms  of  the  higher  harmonies,  i.e.,  q>0  are  exponentially  deereasing  as,  v  > 
kori  and  the  Flankel  funetion  derivative  blows  up,  as  well  as  (-|wi|  |vr|+korisin  \wr\  sinh  |>Vi| 
<0)  and  gives  more  attenuation  with  inereasing  Vrand  fixed  |vi|  as  Vi>kori.  For  the  set  q= 
-1,  the  total  eomplex  phase  is  rewritten  as 
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Figure  6  Contour  deformation  for  shadow  region 
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-jK¥{ y^p  ) «  -|  \{(p  +  h,- 1) + 7(| I  +  (^  + 1 w, I  +  kj^  cos  )  (28) 

the  first  term  shows  that  the  wave  attenuates  as  it  moves  an  angular  distance  (^zi-l-|>Vr|)  with 
an  attenuation  factor  |vi|.  It  is  clear  that  this  term  is  greater  than  the  angle  traveled  by  the 
same  p  wave  but  for  the  set,  q=o,  and  consequently  corresponds  to  higher  attenuation  for 
this  term  and  hence  q=o  is  the  dominant  term  for  the  source  s=o  (see  Fig.  7). 

To  account  for  the  contributions  of  the  souree  s=l  (nl2<^n),  the  eontour  must  be  elosed 
in  the  upper  half  plane  as  shown  in  Fig.  7,  it  is  required  to  replace  (/)hy  {</)  -2n),  hence  the 
total  phase  is  given  by 

-JKwi-yl)  ~ I K |(^ -\wr\-^)- iVi cos w,  +  j\ p, \{(p-7r- 1 wj)  (29) 

where  the  first  term  of  (29)  shows  that  as  if  the  wave  traveling  the  angular  distance  \  (t)-Ti- 

|wr||  with  an  attenuation 
rate  |vi|.  This  angular 
distanee  is  larger  than  that 
of  the  set  (q=o,  s=o),  i.e., 
more  attenuation  is 
experienced,  but  it  is  only 
eomparable  to  it  if  »  n, 
which  will  interfere  with  it 
giving  ripples  in  the 
element  pattern  gain  near 
While  the  seeond 
term  of  (29)  shows  that 
this  wave  has  negative 
phase  veloeity  with 
negative  angle,  so  it  is  a 
forward  wave. 


EXCITED  EUEMENT 


1 


Fig.  7  Ray  interpretation  in  the  shadow 


Similar  considerations  for  the  set  s  =1,  q  =  -1,  with  the  contour  of  integration  is  closed  in 
the  upper  half  plane  shows  that  a  total  phase  given  by 

-jkM-yp)^  y\i(p-^-\w))-j\  cosw,  (30) 


which  shows,  for  the  same  p, 
the  largest  attenuation,  due  to 
the  longest  angle  of  travel 
compared  to  the  three  other 
sets,  i.e.,  (s=o,  q=0),  (s=o,  q=  - 
1),  (s=l,  q=o).  It  is  to  be  noted 
that  the  set  of  creeping  wave 
poles  eorresponding  to  q  =  o 
requires  special  consideration 
since  the  Airy  approximation 
must  be  used  for  the  Hankel 
functions  in  the  vicinity  of 
v»kori.  In  general,  the  residue 
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Figure  8  Contour  in  the  shadow  boundary 
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at  corresponding  to  a  creeping  wave: 


Re^l 


(<L 


=vmlV{^r^e 


-iVp(<p-niT) 


(31) 


where  ap  is  given  by  (21)  given  before,  For  (^zi-7t:/2)>0  (sufficiently  large),  the  residue 
series  is  rapidly  convergent,  since  contribution  from  successive  values  of  p  decay 
exponentially.  However,  when  {(p- n !  i.Q.  in  the  shadow  boundary,  the 

convergence  slows  down  and  an  alternate  technique  must  be  used. 

C. Shadow  Boundary 

Near  the  shadow  boundary  the  harmonic  series  is  poorly  convergent,  and  the  techniques 
used  in  the  lit  region  cannot  be  extended  since  the  Debye  approximation  to  the  Hankel 
functions  is  not  valid  for  v  ~  kori.  This  difficulty  has  already  been  resolved  for  smooth 
cylinders  with  surface  admittance  [1].  We  can  ascertain  that  the  integrand  decays 
exponentially  away  from  v/kori  =  1  along  the  path  F.  Thus,  the  integral  along  F  in  may 
be  approximated,  and  from  the  large  argument  approximation  for  the  Airy  function  [1] 
we  can  ascertain  that  the  integrand  decays  exponentially  away  from  v/kori~/l  along  the 
path  F  (see  Fig.  8).  Thus  the  integral  along  F  in  may  be  approximated  by 

n«>)lr  =  <Ki-,)V‘(Kr„OI  (32) 

where  k  =  n I  T)\ a) 


andg(^,«)=[  e  re  + aAi{2^'^  re  (33) 

vooe  ^ 

The  shadow  boundary  transition  function  g{g,a)  is  a  generalization  of  the  Fock  function. 


V.  Example: 

Single  Mode  Aperture  Field  Approximation 

For  illustration  purposes,  a  TEM  parallel  plate  waveguide-fed  array  is  considered.  The 
aperture  field  is  assumed  to  be  uniform  across  the  aperture.  Matching  the  array  under  the 
equi-phase  condition  v  =  0  (broad-side  scan),  was  selected  for  all  calculations.  The 
numerical  evaluation  of  the  E-plane  element  patterns  was  performed  with  the  help  of 
harmonic  series  for  a  number  of  dielectric  sleeve  thickness  with  8r  =  2.56,  koro  =50  and 
a/?.  =  0.5.  Eigs.  (9),(  10)  show  main  features  of  the  E-plane  element  voltage 
gainpatterns,  normalized  to  the  unit  cell  gain  y2m  /  X  . 

The  presence  of  a  dielectric  cover  on  the  cylindrical  array  is  seen  to  introduce  a 
pronounced  dip,  and  pattern  ripple  near  broadside,  which  tend  to  be  larger  than  those  of 
uncovered  circular  arrays  with  a/A,  >o.  The  location  of  the  pattern  dip  is  essentially 
identical  to  that  observed  in  corresponding  planar  arrays  with  a  dielectric  cover  of  the 
same  thickness.  It  is  observed  that  radomes  of  increased  thickness,  the  dip  moves 
towards  broadside  and  becomes  more  pronounced.  More  ripples  between  broadside  and 
the  dip  are  observed  for  thinner  covers,  while  in  the  shadow  region  almost  a  constant 
slope  (on  a  dB  scale)  is  found.  Characterization  of  these  features  is  obtained  by  the 
asymptotic  valuation  of  the  Integral  representation.  The  results  of  this  calculation  are 
plotted  in  Eigs.  (9)  and  (10)  for  different  radome  thickness  with  koro  =  50. 


In  the  lit  region  0<^zi<90°,  the  total  pattern  is  the  result  of  the  superposition  of  the  space 
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wave  (saddle  point  contributions)  and  contributions  from  fast  wave  poles.  These  poles 
were  considered  in  pairs  {±v^p).  In  numerical  evaluation,  their  number  was  varied  to 

include  just  a  surface  wave  contribution  p=o,  or  in  addition  up  to  four  creeping  waves  (p 
=  0,  1,  2,  3,  4),  the  ripple  in  the  main  beam  in  the  lit  region  is  found  due  to  the 
interference  between  the  space  wave  and  the  fast  set  of  poles  However,  it  was 

determined  that,  only  the  surface  wave  pole  (±Vo*)  contribution  has  a  significant  effect 
in  the  lit  region.  An  approximate  prediction  of  the  location  of  the  dip  can  be  made  by  the 
planar  formula  sin^  =  |v°  /  kfy  -  A,  I  a\,  where  p°  /  is  the  surface  wave  number.  The 
location  of  the  planar  array  nulls  is  indicated  by  arrows  in  Figs.  (9),  and  (10).  The 
pattern  in  the  shadow  region  was  computed  by  the  residue  series  via  only  the  slow  wave 
poles  ( for  (j)>  90°),  where  p  was  varied  from  o  to  3. 


SPACE  WAVE 


angle  (decrees) 

Figure  9  Element  Gain  Pattern. 


Deep  in  the  shadow  region,  the  pattern  is  due  to  the  wave  with  the  lowest  attenuation 
constant,  which  in  the  examples  given  is  the  first  slow  creeping  wave  (p  =  1,  q  =  0,  s  = 
0).  For  ^=Ti,  interference  between  the  two  lowest  creeping  waves  (s=+l,  p=l,  q=0)  and 
(s  =0,  p  =1,  q  =0)  for  yields  the  well-known  back-lobe  ripple  of  the  pattern.  A  detailed 
break  down  of  the  various  wave  contributions  to  the  element  pattern  is  discussed  in  Fig. 
(9).  It  is  seen  that  the  dip  in  the  space  wave  appears  near  44.6°.  For  0<^I5°  and 
50°<^75°,  the  influence  from  the  captured  fast  creeping  wave  poles  and  the  surface 
waves  ((Fj',-Fj')  as  considered  as  their  respective  residue  contributions.  In  the  region 
15°  «^50°,  the  pole  transition  effects  were  taken  into  account. 

As  seen  from  Fig.  9  the  general  outline  of  the  pattern  in  the  lit  forward  region  is 
determined  by  the  space  wave.  By  comparing  the  results  of  harmonic  series  with  that  of 
superposition  of  the  space  wave  and  surface  wave  residue  contribution,  one  concludes  as 
mentioned  that,  the  ripple  in  the  main  beam  may  be  directly  attributed  to  interference 
between  the  space  wave  and  the  fast  surface  wave  pole.  In  addition,  to  obtain  a  good 
accuracy,  one  clockwise  and  one  counterclockwise  fast  creeping  wave  should  be  taken 
into  account.  Again,  the  pattern  in  the  shadow  region  (see  Figs.  9,  10)  was  computed 
sing  residue  series  of  only  slow  wave  poles  (  for  Till«j)<ii),  whose  number  was  also 
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varied  from  p=  o  to  3.  This  produced  sufficient  accuracy  obtained  even  close  to  the 
penumbra. 


Deep  in  the  shadow  region,  the  slope  (in  dB)  of  the  element  pattern  is  essentially  that  of 
the  single  lowest  creeping  wave  pole  contribution.  The  pattern  in  the  shadow  boundary 
transition  75°  <^<110°  was  determined  by  a  numerical  integration  of  the  transition 
function.  Fig.  9  shows  the  element  pattern  for  koro  =50,  element  spacing  a/A,=0.5  Sr=2.56 
and  dielectric  thickness  t/k=0.2,  since  dip  appears  now  at  an  angle  =38°  and  the  pattern  is 
correspondingly  narrower  than  for  t/A,=0.15.  The  dip  level  however,  is  more  pronounced 
than  for  t/X  =0.15.  Also,  the  ripple  frequency  has  increased  as  the  wave  gets  slower. 
Main  factors  in  the  location  of  the  dip  is  the  element  spacing  a/k,  and  the  surface  wave 
location,  as  example,  the  dip  will  appear  more  closer  to  broadside  for  larger  a/A  as 
shown  in  Fig.  (10)  for  a  dielectric  sleeve  of  t/A,=0.15  with  a/A,=0.525.  In  this  case,  the  dip 
occurs  at  an  angle  38°,  compared  to  44.6°  shown  before  for  t/A,=0.15  but  a/A,=0.5. 


For  a  cylinder  of  larger 
radius,  the  above- 
mentioned  effects  are 
modified.  For  the  case  of 
koro  =105  and  a/A,=0.525, 
the  element  pattern  in  the 
lit  region,  is  dominated  by 
the  space  wave,  while 
ripple  amplitude  decreases 
and  its  frequency 
increases  as  for  the  case  of 
cylindrical  arrays  without 
dielectric  radome  [1,2]. 


VI- Conclusions: 

Guided  wave  effects  are  profoundly  influenced  by  curvature,  and  are  clearly 
displayed  in  the  element  pattern  in  mutually  coupled  array.  Element  pattern  serves  as  a 
sensitive  indicator  of  the  degree  of  suitability  of  the  phased  array  radome  combination. 
A  dielectric  sheet  in  contact  with  the  radiating  structure  of  the  phased  array  gives  rise  to 
scan  volume  limited  by  blind  spots  arising  when  phase  synchronization  of  one  of  the 
Floquet  modes  with  guided  surface  wave—  due  to  the  dielectric  cover. 
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